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ABSTRACT 

We examine the effects of a hydrodynamical wind on advection dominated accretion flows 
with thermal conduction in the presence of a toroidal magnetic field under a self-similar treat- 
ment. The disk gas is assumed to be isothermal. For a steady state structure of such accretion 
flows a set of self similar solutions are presented. The mass-accretion rate M decreases with 
radius r as M ex r s+ 2 , where s is an arbitrary constant. We show that existence of wind will 
lead to enhance the accretion velocity. Cooling effects of outflows or winds are noticeable 
and should be taken into account for calculating luminosity and effective temperature of op- 
tically thin and thick ADAFs. Increasing the effect of wind, decreases the disk's temperature, 
because of energy flux which is taken away by winds. We will see that for a given set of input 
parameters, the solution reaches to a non-rotating limit at a specific value of <fi s . The values of 
this limitation on qb s will increase by adding s, wind parameter. In fact, the higher values of s 
will allow the disk to have physical rotating solution for larger <p s . 
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1 INTRODUCTION 

Accretions onto compact objects are the most energetic process in 
the universe. It is believed that many astrophysical objects are pow- 
ered by mass accretion on to black holes. The standard geomet- 
rically thin, optically thick accretion disk model can successfully 
explain most of the observational features in active galactic nuclei 
(AGN) and X-ray binaries (Shakura & Sanyav 1973). 

Accretion disks are an important ingredient in our current un- 
derstanding of many astrophysical systems on all scales. Examples 
of presumed disks accretors include young stars, compact objects 
in close binary system, active galactic nuclei (AGN) and quasars 
(QSOs). There is also evidence that the process of mass accre- 
tion via a disk is often and perhaps always associated with mass 
loss from the disk in the form of a wind or a jet. Mass loss ap- 
pears to be a common phenomenon among astrophysical accretion 
disk systems. These mass-loss mechanisms are observed in micro- 
quasars, Young stellar objet and even from brown dwarfs ( Ferrari 
1998,Bally,Reipurth & Davis 2007; Whelan et al 2005). An outflow 
emanating from an accretion disk can act as a sink for mass, angu- 
lar momentum and energy and can therefore alter the dissipation 
rates and effective temperature across the disk(Knigge 1999). The 
accretion flows lose their mass by the winds as they flow on to the 
central object. As a result of mass loss,the accretion rate, M, is no 
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longer constant in radius r. It is often expressed as M oc r s with s 
being a constant of order unity (Blandford & Begelman 1999). 

Astrophysical jets and outflows emanating from accretion 
disks have been extensively investigated by many researcher. Vari- 
ous driving sources are proposed, including thermal, radiative and 
magnetic ones. Traditionally the name of wind depends on its driv- 
ing force. In this study we will follow the hydrodynamical (ther- 
mal) wind which has been discussed by many authors (e.g. Meier 
1979, 1982, Fukue 1989, Takahara et. al 1989). 

Accretion disk models have been extensively studied during 
the past three decades (see Kato et al 1998 for a review). Besides 
introducing the traditional standard disk by shakura and sunayev 
(1973), there are new-type disks, such as advection-dominated 
accretion flows (ADAF) or radiative inefficient accretion flows 
(RIAF) for a very small mass-accretion disks (Narayan & Yi 
1994) and supercritical accretion disk or so-called slim disks, for 
a very large mass-accretion rate (Abramowicz et al 1998). ADAFs 
with winds or outflows have been studied extensively during re- 
cent years, but thermal conduction has been neglected in all these 
ADAFs solutions with winds and toroidal magnetic field. Since the 
advection-dominated disks have high temperature, the internal en- 
ergy per particle is high. This is one of the reasons why advective 
cooling overcomes radiative cooling. For the same reason, turbulent 
heat transport by conduction in radial direction is non-negligible 
in the heat balance equation. Shadmehri (2008), Tanaka & Menou 
(2006) have studied the effect of hot accretion flow with thermal 
conduction. Shadmehri (2008) has shown that thermal conduction 
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opposes the rotational velocity, but increases the temperature. In ad- 
verted dominated inflow-outflow solutions (ADIOS) it is assumed 
that the mass flow rate has a power low dependence on radius with 
the power low index s, treated as a parameter ( e.g. Blandford & 
Begelman 1999). Kitabatake , Fukue & Matsumato (2002) studied 
supercritical accretion disk with winds, though angular momentum 
loss of the disks, because of the winds, has been neglected. In this 
paper, we present self-similar solutions for ADAFs with thermal 
conduction, wind and a toroidal magnetic field, in the basis of so- 
lutions presented by Akiziki & Fukue (2006). In $2, we present 
assumptions and basic equations. Self-similar solutions are inves- 
tigated in section 3. The aim of this work is to consider possibility 
that winds in the presence of thermal conduction and a toroidal 
magnetic field which has been largely neglected before, could af- 
fect the global properties of hot accretion flows substantially. In $4 
we show the results and give discussions of the results. 



2 THE BASIC EQUATIONS 

We investigate the effect of mass outflow by the wind and mass 
accretion rate by viscosity simultaneously. We consider an accre- 
tion disk that is stationary and axi-symmetric {-Br — 4t — 0) and 
geometrically thin j- < 1. In cylindrical coordinates (r, if, z), we 
vertically integrate the flow equations, also we suppose that all flow 
variables are only a function of r. We ignore the relativistic effect 
and we use Newtonian gravity. We adopt a-prescription for vis- 
cosity of rotating gas. For magnetic field we considered a toroidal 
configuration. 

The equation of continuity gives 



dr {VljVr)+ 2K dr 



= 



(1) 



where v r is the accretion velocity (v r < 0) and £ = 2pH is the 
surface density at a cylindrical radius r. The mass loss rate by out- 
flow/wind is represented by M w . So 

M w (r) — / 4irr'rri w (r')dr' 

where rn w (r) is mass loss rate per unit area from each disk face. 
On the other hand, we can rewrite the continuity equation: 



-^-(rXv r ) = 2pH 
r or 



(2) 



where p the mass loss rate per unit volume and H is the disk half- 
thickness. 

The equation of motion in the radial direction is: 
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where v^ , C 3 and Ca are the rotational velocity of gas disk, sound 
and Alfven velocities of the fluid respectively. Sound speed is de- 
fined as C s = ^22± 5 p gas being the gas pressure and Alfven ve- 
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locity is defined as C\ = ^— 



where p mag being the 



magnetic pressure. 

The integrated angular momentum equation over z gives: 
(e.g., Shadmehri 2008) 
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where the last term on the right side represents angular momen- 
tum carried by the outflowing material. Here, I — corresponds 



to a non-rotating wind and I — 1 to outflowing material that car- 
ries away the specific angular momentum (Knigge 1999). Also v is 
kinematic viscosity coefficient and we assume: 



v = aC s H 



(5) 



where a is a constant less than unity. By integrating over z of the 
hydrostatic balance ,we have: 
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where /3 = J"° a = ^(tt 1 ) 2 which shows the important of mag- 
netic field pressure compared to gas pressure. We will show the 
dynamical properties of the disk for different values of /3. Now we 
can write the energy equation considering cooling and heating pro- 
cesses in an ADAF. We assume the generated energy due to viscous 
dissipation and the heat conduction into the volume are balanced by 
the advection cooling and energy loss of outflow. Thus, 
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where the second term on right hand side represents energy transfer 
due to the thermal conduction and F s — 5$ s pCg is the saturated 
conduction flux (Cowie & Makee 1977). Dimensionless coefficient 
$ s is less than unity. Also, the last term on right hand side of energy 
equation is the energy loss due to wind or outflow (Knigge 1999). 
Depending on the energy loss mechanisms, dimensionless parame- 
ter rj may change. In our case we consider it as a free parameter in 
our models so that larger r\ corresponds to more energy extraction 
from the disk because of outflows (Knigge 1999). Finally since we 
consider the toroidal component for the global magnetic field of 
central stars, the induction equation with field scape can be written 
as: 

d 
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where B v is the field scaping/creating rate due to magnetic insta- 
bility or dynamo effect. 



3 SELF-SIMILAR SOLUTIONS 

Self-similar solutions can not be able to describe the global be- 
havior of the solutions, because in this method there are not any 
boundary conditions which have been taken into account. However 
as long as we are not interested in the solutions near the boundaries, 
such solution describe correctly, true and useful asymptotically be- 
havior of the flow in the intermediate areas. 

We assume that the physical properties are self-similar in the 
radial direction. In the self-similar model the velocities are assumed 
to be expressed as follows: 
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and constants C\,Ci and C 3 will be determined later. From the 
hydrostatic equation, we obtain the disk half-thickness H as: 



H 



— = VC 3 (l + /3)=tan< 



(14) 



Hence, a hot accretion flow with winds also has a conical surface, 
whose opening (half-thickness) angle is a. 

We assume the surface density E is in the form of: 



S r s 



(15) 



Note that the value of s should be determined iteratively for 
consistency. Then we assume that the power law index of the den- 
sity p in the radial direction is constant regardless of z. Hence we 
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It should be noted that, for a self-similar disk without any wind 
mass loss, the suffix sis s — — |. 

By substituting the above self-similar solutions into the dy- 
namical equations of the system, we obtain the following dimen- 
sionless equations, to be solved for Ci, C'2 and C 3 : 
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where m = 



Mqw 



is the nondimensional mass accretion rate. 
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After algebraic manipulations, we obtain a fourth order algebraic 
equation for C\ : 

D 2 Ct + 2DBCI + (B 2 + 2D(E - lj)Cf + 

(2B(E - 1) - A 2 )Ci + (E - l) 2 = 

Where 
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Figure 1. Numerical coefficient Ci s as a function of advection parameter / 
for several values of 0, magnetic field strength. All of this figures was set 
up for s = -0.3, a = 0.01, <j> = 0.001 and I = r\ = 1. 



Abbassi et al. (2008) have solved these equations when s — 
— i, because they did not have any wind or mass-loss in their 
model. This algebraic equation shows that the variable C\ which 
determines the behavior of radial velocity depends only on the a, 
$ s , /3, /, s and r\. Other flows quantities such as C2 and C 3 can be 
obtained easily from C\ : 
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We can solve these simple equations numerically and clearly 
just physical solutions can be interpreted. Without wind, thermal 
conduction and magnetic field, s — I — rj = <p s = /3 = 0, 
the equations and their similarity solutions reduce to the standard 
ADAFs solution ( Narayan & Yi 1994). Also without wind they 
reduce to Abbassi et al. (2008). Now we can analysis the behavior 
of solutions. 



4 RESULTS 

Now we can analysis the behavior of solutions in the presence of 
outflow, thermal conduction and magnetic field. Our primary goal 
is to investigate these effects via parameters s, I, rj, (f> s and j3. 

This Algebraic equation shows that the variable C\ which de- 
termines the behavior of radial velocity depends only on the a, <f>, 
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Figure 2. Numerical coefficient C; s as a function of advection parameter 
/ for several values of <f> s , thermal conduction parameter. All of this figures 
was set up for s = —0.3, a = 0.01, /3 = 0.1 and I = 77 = 1. 



Figure 3. Numerical coefficient d s as a function of advection parameter 
/ for several values of s, wind parameter. All of this figures was set up for 
4> s = 0.001, a = 0.01, 13 = 0.1 and / = r; = 1. 



/3, s and /. Other flows quantities such as C2 and C3 can be ob- 
tained easily from C\ via the above equations. The parameters of 
the model are the ratio of specific heats 7, the standard viscose pa- 
rameter a, the energy-advection parameter /, degree of magnetic 
pressure to gas pressure /3 and s which determines the outflow from 
the disks. 

Figure 1. shows the coefficients d& in terms of advection pa- 
rameter / for different values of /3. In the upper panel we have 
shown C\ as a function of advection parameter. It represents the 
behavior of radial flows of accretion materials. Although in ADAFs 
the radial velocity is generally sub-Keplerian, it becomes larger by 
increasing /. By adding /3, which indicates the role of magnetic 
field on the dynamics of accretion disks, we will see that the ra- 
dial flow will increases. On the other hand the radial flow increases 
when the toroidal magnetic field becomes large which agrees with 
Abbassi et al. (2008)'s results. This is because, the magnetic tension 
term dominates the magnetic pressure term in the radial momentum 
equation, which assists the radial velocity of accretion flows. The 
middle panel, shows the coefficient C2 which is the ratio of the ro- 
tation velocity to the Keplerian velocity as a function of advection 
parameter for different values of /?. We can see that by adding the 
magnetic field influences (adding /3), the rotation velocity will de- 
creases. This is because of that the disk should rotate faster than the 
case without the magnetic field which results the magnetic tension. 
In the lower panel, we have plotted C3, the coefficient of squared 
sound speed as a function of /. As the advection parameters be- 
comes large, the sound speed and therefore the disk's thickness be- 
comes large. Sound speed also depends on the magnetic field influ- 
ences. By adding magnetic field influence we will see that sound 
speed increases as well as vertical thickness of accretion flows. 



In Fig. 2 we investigate the role of saturated thermal conduc- 
tion in radial, toroidal and sound speed of accretion flows. Increas- 
ing the thermal conduction coefficient <t> s , will decrease the radial 
velocity. It will have a large effect on the rotational velocity of ac- 
cretion flows (Middle panel). It will decrease the sound speed and 
therefore vertical thickness of the disks too (lower panel). 

To show the behavior of the solution respect to the wind influ- 
ences we have plotted disk's physical quantities for different values 
of s and / in Fig.3 and Fig.4. The value of s measures the strength of 
outflows, and larger values of s denotes strong outflows. Each curve 
is labeled with its corresponding s. We can see that ADAFs with 
wind rotates more quickly than those without wind, radial flows 
of the accretion materials also increases for larger s. Strong-wind 
models have a lower vertical thickness compare to the weak-wind 
model because by adding s we can see that C3 decreases which 
means that the sound speed and vertical thickness decreases as well. 

There are some limitations for having hot disks with sufficient 
low /. Middle panel in Fig.3 shows that there is no physical so- 
lution for ADAFs with / < 0.4 for the case of no-wind. But by 
adding s we can see that this limitation can be canceled. The out- 
flow and wind play as a cooling agent while thermal conduction 
provides an extra heating source so these two factors affects the 
dynamical behavior of the flows. While effects of wind and ther- 
mal conduction on the profiles of sound speed and vertical thick- 
ness of the disks is similar, their effects on the profiles of radial and 
rotational velocities oppose each other. Both of them leads to de- 
crease the sound speed and therefore vertical thickness of the disks 
decreases. They lead to enhanced radial and rotational velocities. 

In the lower panel of Fig.3 and Fig.4 we have plotted C3 for a 
sequence of /. They show two cases, / = r\ = and I — r\ — 1. In 
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Figure 4. Numerical coefficient d s as a function of advection parameter 
/ for several values of s, wind parameter. All of this figures was set up for 
4> s = 0.001, a = 0.01, P = 0.1 and I = 0, ry = 0. 



Figure 5. Numerical coefficient d s as a function of thermal conduction 
parameter <f> s for several values of 0, magnetic field parameter. All of this 
figures was set up for / = 0.7, a = 0.001, s = —0.3 and I = 1, rj = 1. 



the case of / = r\ = 0, variations of s do not affect on Cz profiles 
( and vertical thickness of the disks); while in the case I — r) — 
1 it decreases when we increase the wind effect. Higher r\ means 
more energy flux is taking away by winds and I will account for 
the angular momentum carry out by the wind. So solutions with 
no-rotating wind I — and no energetic wind r\ = was plotted in 
Fig.4. 

In Fig. 5, Fig.6 and Fig.7 we have plotted the coefficient C\s 
with thermal conduction parameter cf> s for different values of j3, f 
and s. Tanaka & Menou (2006) have shown that for a very small 
(f> s their solutions coincide the original ADAFs solution Narayan 
& Yi (1994); But by adding saturated conduction parameter, <j> 3 , 
the effect of thermal conduction can be better seen when we ap- 
proach 4> s ~ 0.001 — 0.01. So we have plotted our solutions in 
this range. We can see that for a given set of input parameters, the 
solution reaches to a non-rotating limit at a specific value of <p 3 . 
We can not extend the solution for values of this limit for <fi 3 . Be- 
cause Eq.31 gives a negative value for Cf which is clearly unphys- 
ical. The values of this limitation on (j> s will increase by adding s, 
wind parameter. In fact, the higher values of s will allow the disk 
to have physical rotating solution for larger <fr s , Fig.7 middle panel. 
In Fig.5 we can see that magnetic field has an important role on the 
radial velocity and the sound speed of the flow but it can not modify 
the rotation velocity which are compatible with Akizuki & Fukue 
(2006), Fukue (2004) and Abbassi et al (2008). The sound speed 
and radial infall velocity both increase with magnitude of thermal 
conduction. 

As the level of thermal conduction is increased, more and 
more heat flows from hot inner regions, resulting in a local increase 
of the gas temperature relative to original ADAFs solutions. Simul- 



taneity, the gas adjust its angular velocity and increases its inflow 
speed to conserve its momentum balance. We have found that the 
level of advection, /, will modified the dynamical structure of disk 
significantly. The breakdown the solutions, when Cl — > 0, occurs 
at lower values of /. Clearly Thermal conduction can significantly 
affect the structure and the properties of the hot accretion flows. To 
better understand the role of thermal conduction it is better to study 
it in a more realistic two-dimensions investigations, without height 
integrations. 

Our results that shows out flow (large s) in ADAFs will mod- 
ify the dynamical quantities of accretion flows significantly. Re- 
cently Kawabata & Mineshige (2009) have shown that thermally 
driven wind with a power law formula is expected for such hot 
ADAFs disks. They have also shown that the ratio of the outflow 
to the accretion rate is approximately unity for a ^ 0.1. So our 
investigation should have significant role to understand the physics 
of radiation inefficient accretion flows RIAFs or ADAFs. 

However such kind of solutions enable us to investigate easily 
all possible cases for hot accretion flows with thermal conduction, 
wind, B-field and any combinations of these physical mechanisms. 
This simple toy model will prepare a frame work to compare these 
effects. 



5 GENERAL PROPERTIES OF ACCRETION FLOWS 

In the pervious sections we have introduced solutions of advection 
type accretion flow with winds in the presence of thermal conduc- 
tion bathed in a toroidal magnetic field. In this section we will dis- 
cuss the general properties of accretion flow and we will investigate 
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Figure 6. Numerical coefficient ds as a function of thermal conduction 
parameter 4>s for several values of /, advection parameter. All of this figures 
was set up for a = 0.01, /3 = 0.1, s = — O.s and I = 1, 77 = 1. 



Figure 7. Numerical coefficient d s as a function of thermal conduction 
parameter (f> s for several values of s, wind parameter. All of this figures 
was set up for a = 0.01, /3 = 0.1,/ = 0.7 and I = 1, tj = 1. 



possible effects of wind, thermal conduction and magnetic field on 
the physical quantities of accretion flows. 

Using the self-similar solution we will estimate the mass- 
accretion rate as: 



M = -27rrEw r = 27rE ciaVGMr s+1/2 
r ns+i/2 



= M out { — y 



f*out 



(33) 



where r ou t is the disks outer radius and M ou t is mass-accretion 
rate there. In the case of accretion disk with no wind, s = —1/2 
the accretion rate is independent of the radius, while for those with 
wind, s > —1/2, the accretion rate decreases with radius as we 
expect. Because winds start from various radii, the mass loss rate 
is not constant but depends on the radius. As a result, some parts of 
accretion materials are not concentrated at the center, but are dilated 
over a wide space. According to wind mass loss, the accretion rate 
decreases with radius as we expected. 

With a simple calculation we will see that a significant amount 
of accretion materials flung away via the wind and mass loss. Only 
1 — 10% of M ou t is ultimately accreted onto the central accretor. 
Using the above expression for accretion rate, we have: 



M in 



1 r in ss+1/2 



M ou t r out 

if we estimate n n ~ W~ 3 r ut we finally have for s — 1/2: 

Min _ 10 -3 

M out 

We can show the radial behavior of temperature of the present 
self-similar ADAFs disks with toroidal magnetic field and outflow. 



These advection-dominated accretion flows occurs in two regimes 
depending on their mass accretion rate and optical depth. Opti- 
cal depth of accretion flows are highly depends on their accretion 
rate. In a high-mass accretion rate, the optical depth becomes very 
high and the radiation generated by accretion flows can be trapped 
within the disk. This type of accretion disks is known as optically- 
thick or slim disks. In the limit of low-mass accretion rates, the disk 
becomes optically thin. In this case, the cooling time of accretion 
flows is longer than the accretion time scale. The energy generated 
by accretion flows therefore remains mostly in the disks and the 
disks can not radiate their energy efficiently. We will examine the 
influence of the wind and thermal conduction on the radial appear- 
ances of the temperature in these two cases. In the optically thin 
case, where the gas pressure is dominated, we can adopt the ideal 
gas law to estimate the effective temperature as (Akizuki & Fukue 
2006): 

GM 






c s 



(34) 



where T is the gas temperature, R gas constant and ft the mean 
molecular weight. If we use (GM — r g c), where r g and c are the 
Schwartzshild radius and light speed respectively, the temperature 
gradient is expected as: 



T = c 3 



c 2 /2 r g 
R/Jx r 



(35) 



which means that T oc — . This has similar form ( radial depen- 
dency) with non magnetic case, but coefficient C3 implicitly de- 
pends on magnetic field j3, outflow effect s, advection parameter / 
and the effect of thermal conduction cj>. For the case of strong wind 



Cz decreases, Fig. 3, so the surface temperature and vertical thick- 
ness of the disk will decrease. Thermal conduction has the same 
effect, Fig.2, while the magnetic field will increase the surface tem- 
perature and vertical thickness of the disk. 

It should be emphasized that the radiative appearance of op- 
tically thin disks such as L v , can not be calculated easily. In this 
case it should demonstrate the importance of advective cooling in 
optically-thin disks. In order to clarify this problem we should show 
that cooling and heating how change with r, E and M. Radiative 
cooling generally has very complicated parameter dependencies. In 
the optically thin cases, emission from the gas is not a blackbody 
continuum. Bremsstrahlung cooling by non-relativistic or relativis- 
tic, Synchrotron and Compton cooling may have possible role to 
reproduce emission spectra. 

In the optically thick case, where the radiation pressure dom- 
inated, sound speed is related to radiation pressure. We can write 
the average flux F as: 



F = aT c 4 = g^n 9as = -cE 



C3 



1 + 13 



GMr s 



(36) 



where II sas = T,c? s is the height-integrated gas pressure, T c disk 
central temperature and a the Stefan-Boltzman constant. In the 
moderately to the strongly magnetized /3 ^ 1 cases the magnetic 
pressure is comparable with gas pressure and it should take into 
account: 

n = u gas + n mag = (l + /3)ii gas 

so the calculation should be modified by multiply to 1 + /3. For 
optically thick case the optical thickness of the disk in the vertical 
direction is: 

1 v 1 v • 
T — -KE = -KL r , 

where k is the electron-scattering opacity. So we can calculate the 
effective flux and effective temperature of the disk surface as: 



oT, 



oTt 3c / c 3 GM 
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4k V 1 + P r 2 4 V 1 + /3 47rr 2 
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\l/4 r -l/2 



(37) 



(38) 



where Le = 47rc^— is the Eddington Luminosity of the central 
object. If we integrate these equations radially we have the disk 
luminosity as: 



L d 



C3 r -1 Tout 

-——L B In , 

1 + p Tin 



(39) 



As this equation shows, the optically thick disk's luminosity is 
affected by magnetic field explicitly, /?, but it would be affected by 
outflow, thermal conduction and viscosity through the C3 implic- 
itly. 

It should be emphasized that the Luminosity and effective 
temperature of the disk, Ldisk and T e ff, are not affected by the 
mass loss through outflow (there are not any s dependencies). But 
wind would affect on the radiative appearance of the disk through 
the C3 in these formulas, implicitly. The average flux decreases all 
over the disk when we have mass loss outflow compared with the 
case of no mass loss. The surface density and there fore optical 
depth decrease for the mass loss case. So we can see the deep in- 
side of the disk. 



6 SUMMERY AND CONCLUSION 

In this paper we have studied accretion disks around a black hole 
in an advection-dominated regime in the presence of a toroidal 
magnetic field. Thermal conduction and wind as a energy sources 
and angular momentum transport were adopted. We have presented 
the results of self-similar solutions to show the effects of various 
physics behavior on the dynamical quantities of accretion flows. 
We adopt the solution presented by Akizuki & Fukue (2006), Fukue 
(2004) and Abbassi et al (2008) to present dynamical behavior 
of the advection dominated accretion flows. Some approximations 
were made in order to simplify the main equations. We assume an 
axially symmetric, static disk with a-prescription of viscosity. We 
ignored the relativistic effects and the self-gravity of the disks. Con- 
sidering weakly collisional nature of hot accretion flow ( Tanaka 
& Menou 2006, Abbassi et al. 2008), a saturated form of thermal 
conduction was adopted as a possible mechanism for energy trans- 
portation was adopted. We have accounted for this possibility by 
allowing the saturated thermal conduction constant, <f> a , to vary in 
our solutions. 

Theoretical arguments and observations suggest the mass loss 
via wind in RIAFs. By assuming that the flow has self-similarity 
structure in the radial direction, a power-law wind is adopted in our 
model. Using some assumptions, we made a simplified toy model 
which included the effect of thermal conduction, wind and B -field. 
In this toy model, we can easily investigate its possible combine 
effects on the dynamical quantities of the fluid. 

We have shown that strong wind could have lower temperature 
which are satisfied with the results were presented by Kawabata & 
Mineshige 2009, which could make significant differences in the 
observed flux compared to standard ADAFs. The most important 
finding of our simple self-similar solutions is that the accreting flow 
is strongly affected not only by mass loss but also by energy loss 
by the wind. 

There are some limitations in our solutions. One of them is 
that the magnetic field also has an important role to produce wind 
(magnetically driven wind). To having X-wind, pure toroidal B- 
field in not enough and the disk should have a z-component B- 
field. So our case is not a good model for magnetically driven disks. 
The other limitation of our solutions is the anisotropic character of 
conduction in the presence of magnetic field. Balbus (2001) has 
argued that the dynamical structure the hot flows could be affected 
by anisotropic character of thermal conduction in the presence of 
B-field. 

Although our preliminarily self-similar solutions are too sim- 
plified, but clearly improve our understanding of the physics of 
advection-dominated accretion flows around a black hole. For hav- 
ing a realistic picture of accretion flow a global solution needed 
rather than the self-similar solution. 
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